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Abstract 


Studies of meson spectroscopy have often employed a non- relativistic Coulomb plus 
Linear Confining potential in position space. However because the quarks in mesons move 
at an appreciable fraction of the speed of light, it is necessary to use a relativistic treatment 
of the bound state problem. Such a treatment is most easily carried out in momentum 
space. However the position space Linear and Coulomb potentials lead to singular kernels 
in momentum space. Using a subtraction procedure we show how to remove these sin- 
gularities exactly and thereby solve the Schrodinger equation in momentum space for all 
partial waves. Furthermore, we generalize the Linear and Coulomb potentials to relativis- 
tic kernels in 4-dimensional momentum space. Again we use a subtraction procedure to 
remove the relativistic singularities exactly for all partial waves. This enables us to solve 
3-dimensional reductions of the Bethe-Salpeter equation. We solve six such equations for 
Coulomb plus Confining interactions for all partial waves. 


PACS numbers: ll.10.St, ll.10.Qr, 14.80.Dq 
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1. INTRODUCTION 


Meson spectroscopy 1 has been one of the most interesting and fundamental subjects 
in elementary particle physics for the last two decades. It has provided one of the basic 
testing grounds for our understanding of both the symmetries and the dynamics of the 
strong interaction between quarks, mediated by gluons. Future studies are also of great 
interest, particularly as they may provide evidence of constituent glue. Given the important 
role of meson spectroscopy it is vital that our theoretical descriptions of these relativistic 
qq systems be as accurate and consistent as possible. Thus one would ideally like to be able 
to connect the theoretical description of mesons to the fundamental theory of the strong 
interactions, namely Quantum Chromodynamics (QCD). However the non-abelian nature 
of QCD leads to strong self-interactions between the gluons resulting in field equations that 
are highly nonlinear and are unable to be solved by standard diagrammatic methods except 
in the perturbative regime. In the region of large distances , lattice gauge calculations 2 , 
which provide the most direct link to QCD, have led to the conclusion that in the static 
quark limit the force between quarks can be very well described with a Linearly rising 
plus Coulomb potential. Nonrelativistic models which use such a potential have been 
very successful in accounting for both the masses and decays of mesons, particularly those 
containing heavy quarks. 

However the pure non-relativistic model calculations have limitations. Firstly for 
systems containing one light quark the use of pure nonrelativistic formalism is obviously 
unjustified. Secondly the nonrelativistic formalism has intrinsic problems such as the 
incorrect dependence of the meson mass on the quark mass, i.e. the mesons with light 
quarks can become heavier than the mesons with heavier quarks 3,4 . Also the nonrelativistic 
Linear potential does not lead to Linear Regge trajectories 3 . None of these problems 
occur in semirelativistic treatments where the relativistic expression for the energy is used. 
Clearly then one must also introduce relativistic effects. Such studies have been made 
and good descriptions of the entire meson family have been obtained 5 . However, if one 
incorporates relativity into a position space calculation then many different relativistic 
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effects must be put in ”by hand” leading to a significant number of adjustable parameters 5 . 
A much more satisfactory approach can be made by doing calculations in momentum space 
where relativistic effects can be handled in a much more economic way. 

Such calculations immediately present two difficulties. Firstly, because one would 
like to retain a manifestly covariant approach it is natural to transform the Linear plus 
Coulomb potentials to momentum space. The problem is however, that both potentials 
lead to singular kernels. Secondly, because many mesons of interest contain quarks of 
comparable mass, one should ideally solve the two-body Bethe-Salpeter (BS) equation 6 
and certainly not consider only the one-body Klein-Gordon or Dirac equations. Although 
the best way to do meson physics in the two-body framework would be to solve the Bethe- 
Salpeter(BS) equation, it is more practical to solve a three-dimensional 7-11 reduction of it. 
However there exist in principle infinitely many possible three-dimensional reductions 7-11 
of the BS equation and generally there is no reason to prefer one reduction to another, 
although in some cases the physical problem itself might suggest a particular reduction 
scheme. Therefore for the general qq problem it is useful to carry out a systematic study 
of the various reductions of the BS equation. 

In this paper we present a complete study of how to solve relativistic two-body bound 
state equations in momentum space with kernels which are a generalization of Coulomb 
plus Linear potentials. A method for treating linear and Coulomb potentials in momentum 
space for the nonrelativistic case was presented by Spence and Vary 16 but their method is 
not easily generalizable for the relativistic case if one retains retardation in the interaction. 
We present a systematic treatment of how to deal with the momentum space singularities 
for both Coulomb and Confining interactions for all partial waves and for both the non- 
relativstic Schrodinger equation and for six different 3-dimensional reductions 7-11 of the 
Bethe-Salpeter equation. The only parameters that our method permits are the quark 
masses and the Coulomb and Confining couplings. Our study is a comprehensive treatment 
of relativistic effects but with a very restrictive parameter set and should thus eventually 
provide a definitive description of the entire meson spectrum. The main purpose of the 
present paper is to present the theoretical subtraction techniques necessary to solve two-body 
relativistic bound state equations in momentum space. 
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2. SCHRODINGER EQUATION AND POTENTIALS IN MOMENTUM SPAC 


The nonrelativistic power law potential in r-space can be written as 


( 2 . 1 ) 

i|— »o v ' 

Where A ^ is the strength of the potential and tj is the screening parameter. The 
index N indicates the type of potential under consideration, i.e N = — 1 corresponds to 
Coulomb potential (A;v = Ac) and N = 1 corresponds to a Linear potential (Ayv = A^,). 
In the present paper we shall be considering only these two types of potential. For the 
bound state problem of two particles with masses mj and m 2 interacting via F^q) the 
Schrodinger equation in momentum space is 

|^(p) + J V" v (qMp')V = BK p) (2.2) 

where p is the reduced mass. The momentum space potential is given by the Fourier 
transform of Eq.(2.1) namely 
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Where q = p' — p. The Schrodinger equation for the I th partialwave is given by 


( 2 . 3 ) 


2 °° 

+ J Vl N {p',p)<i>r,l(p')p' 2 dp' = E n i4> n i{p) 


( 2 . 4 ) 


where p = |p|, n is the principal quantum number and l is the orbital quantum number. 
The partial wave components of the potential is readily obtained as 

1 

T r N / \ n / \ 

lx 


F/V,p) =2* J V N (q )P,(x)d 2 
-*1 

7T r)— *u' dr fN+1 ppt 


where x = cos9 pp > and y is defined as 


V = 


P 2 + P >2 + TJ 2 

2 PP 1 


( 2 . 5 ) 


(2.6) 


5 


The exact lim,_u will be taken shortly. Special cases of interest are the Coulomb case 
( N = —1) and the Linear potential ( N = 1) and they are readily obtained from Eq(2.5) as 


and 


Ac Qijy) 


V, C (p',p) = — Km. 


7 r n—o pp 1 


(2.7) 


v, l (p,p) 


Al ,. d 2 Qi(y ) 
— 11111 7~ 

7 T CT]* pp' 


— lim 

7 r i/— o 


QVy) , ± 

.{pp 1 ) 2 (pp') 


■Q"i(y)\ 


(2.8 a) 
(2.86) 


Here Q;(y) are the Legendre polynomials of the second kind and their first and second 
derivatives are taken with respect to y, i.e. 


«',(») = (2-9a) 

Q",(V) = (2.96) 

We note here that these potentials (at 7 / = 0 ) have singularities when p f — p which 
corresponds to y = 1. In order to see the singularity structure explicitly we rewrite Qt(y) 
in terms of Qu(y) as 


where 


Note also that 


Qi(y ) = Pi(y)Qu(y) - ™i-i(y) 

‘ x 

™i-i(y) = — Pi-m(y)Pm-\(y) 

* ■* TTL 


Qu(y) = 1/2 /n|(y + l)/(y - 1)| = 1/2 In 


V+p) 2 + t, 2 ' 

,(p' ~P ) 2 + r? 2 . 


(2.10a) 


(2.106) 


(2.11a) 


Qo(y) 



-1 


= PP 


+ 


(p' - P ) 2 + y 2 (p' + p ) 2 + y 2 J 


(2.116) 
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and 


1 1 Q'M = r, 2 (p 2 +P°+V 2 ) 


PP 


-1 


+ 


[{p'-pY+v 2 {p'+pY + v 2 J 


( 2 . 11 c) 


In the expression for Qi(y) the only term that is singular (at tj — 0 ) is Qo(y)* Therefore 
the Coulomb potential has a logarithmic singularity from Qo{y) and the Linear potential 
has higher order singularities from Q'o(y) and Q" o(y)* We note that the singularity struc- 
ture of these potentials are the same for all partial waves. 

As mentioned above the potentials we are interested in have singularities at p' = p 
and in the following sections we will show how to take care of these singularities in the 


momentum space Schrodinger equation. There are two useful integrals which will be used 


repeatedly in the following sections. They are 


/ 


<?u(y,q = Q) , / _ ^ 
p' 2 


f[—,Q\{y) + <?\>(y )W = o 
J pp 

0 


( 2 . 12 ) 


(2.13) 


2.1 NON-RELATIVISTIC COULOMB PROBLEM 


In this subsection we will present a subtraction method which will treat the Coulomb 
singularity properly. For the pure Coulomb problem in momentum space the exact analytic 
bound state solutions were found by Fock 12 , but our aim is to solve the Schrodinger 
equation and later relativistic equations for a combined Linear plus Coulomb interaction. 
Thus we need to be able to implement a numerical subtraction procedure in momentum 
space. Apart from the rearrangement of terms this method is identical to the one developed 
in references 13 and 14, but we reproduce it here for completeness. With the potential 
given in Eq.(2.7) and using the expression ( 2 . 10 a) for Qi(y) the Schrodinger equation (2.4), 
with 77 = 0 becomes 
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oo oo 

77-4>nl{p) + — [ Pl{y)^~^-<f>nl{p')p 2 dp' - — [ V>l-l(y)<t>nl(p')p' d P = E n l<t>nl(p) 

2/1 7T p J p ' 7T p J 

0 0 

(2.14) 

Since tu/_i(y) contains no singularity the second integral needs no special treatment. In 
order to remove the singularity arising from the Qo(y) term we subtract and add a term 
from the integrand fo the first integral of Eq(2.13). The added term is propotional to the 
integral of Eq(2.12) and we obtain a singularity free equation 


£ m ?) + *£ 

2/1 7T p 


oo 

J Pi(y) q y 1 [p' 2 Mp')~ 


p 2 4ij>i(p) 

PM 


dp' 


+ ~~ [yP 2 <Mp) 

7T p [ 2 


Ac 
t rp 


OO 

J tn/-l(y)0n«(p')p ,<i P = Enl<t>nl{p ) (2.15) 


Note that at the singular point we have p = p' , y = 1 and Pj( 1) = 1. Therefore the terms 
in the square brackets cancel exactly and removes the singularity arising from (?u(y)- The 
numerical solution of this equation is discussed in the section on Numerical Methods (Sec. 
2.3). 


2.2 NON-RELATIVISTIC CONFINING PROBLEM 


In the case of the Linear potential there are singularities arising from Qu(y)> Q'uiy ) 
and Q" u (y). We are interested in solving the Schrodinger equation in the limit 77 = 0. For 
the sake of clarity we will first consider the / = 0 case. For 1 = 0 the potential is 


K/V.p) = + 


<?'u(yh 

(pp') 2 > 


Therefore the s-wave Schrodinger equation is 


(2.16) 


t-KM + % lim / (^-Q"„(y) + Mp) (2-17) 

2/i 7 rp 2 0 J v pp! 
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Now by adding and subtracting a term in the integral we obtain 


|^<£nu(p) + Kin J {~j^Q"M + Q'M) (^no(p') - 0no(p))dp' 

* 0 

oo 

+ (p) lim f(^Q"„(y) + Q'M) = En »Mp) (218) 

7rp J v—o J pp 1 

o 

From Eq(2.13) we see that the last integral is identically zero. Now we can take the rj = 0 
limit explicitly 15,lc and we finally get 

OO 

^<^nu(p) 4- ^ J <? , o(y)(^no(p) - 0no(p))dp' = £ n O0 n u(p) (2.19) 

0 

In the above equation <?' u (y) has a double pole singularity at p = p' (see Eq2.11b). By 
Taylor series expansion of <t> n u{p ' ) around the point p 1 = p we can see that only a principal 
value singularity is left, which can be treated by conventional means. (see section on nu- 
merical methods and reference 17. Next we consider the case for general l. After removing 
the terms which can be shown to vanish when the T) = 0 limit is taken the form of the 
potential is 


v,Hp'.p) = v , u 5„ 


>(y ) + 


ym \ 

pp 1 ) 2 ) 


+ 


p,[y) 

Substituting this into the Schrodinger equation (2.4) we have 

OO 

Km J \ p M{~iQ n M + Q'M ) 


P'i{y)Qo{y) - w'f-i(y) 


(pp') 2 


( 2 . 20 ) 


4>nl{p ) H ^ 

2/i ?rp z I»- 


+ P'i{y)QM -w'z_i(y) 4>ni{p)dp = E„i<j> n i(p) (2.21) 


In order to remove the singularities now we must perform two separate subtractions. 
The first subtraction is for the singularities coming from Q'u(y) and Q" ^{y) and the second 
subtraction is for Qo{y). For the singularities arising from Q' and Q" , by using Eq.(2.13) 
we can make a subtraction without having to add anything back and for the singularity 
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arising from Qu , by using Eq.(2.12) we subtract and add a term as for the Coulomb case. 
In addition, we have shown previously 15 how to take the explicit tj = 0 limit. Thus our 
singularity free equation, in the exact tj = 0 limit is 


°o 

^nt(p) + ^ fPl(y)Q' o(v) (W) - dp' 


+ 


\l_ 

7rp 2 


OO 


f(f + 1 )p<Mp) 

2 P'i(y) 


dp' 


+ 


A L 1(1 + 1) 


7Tp 


2 2 


-P<t>nl(p)\ 


- — J [ w't-l(y)4>nl(p')dp = E n l4> n l(p) 

*P 2 J 


( 2 . 22 ) 


At the singular point (p 1 = p; y = 1) we have P t (y = 1) = 1 and the bracketed term in the 
first integral vanishes. Therefore, as in the / = 0 case we are left a pricipal value singularity. 
In the second integral at the singular point P'i(y = 1) = /(/ + 1 )/2 and the term in the 
bracket again vanishes and kills the logarithmic singularity arising from Qu(y). Note that 
for / = 0 Eq(2.22) reduces to Eq(2.19). Now we are in a position to solve Eq.(2.15) for the 
pure Coulomb or Eq.(2.22) for the pure Linear potential for all partial waves. It is also 
obvious how to treat the combined Coulomb plus Linear potentials together. 


2.3 NUMERICAL METHODS 


Consider first the Coulomb equation (2.15). An important point to note is that at 
the singular point p' = p the term in the square brackets of the first integral goes to zero 
faster than the logarithmic singularity in Qo(y), and therefore the integrand of the first 
term is identically zero at p' = p. By using Gaussian quadrature one can easily write the 
whole equation (2.15) as a matrix equation with <f n i(p) as the eigenvectors and E„i as the 
eigenvalues. Because the kernel of the first integral is zero when p' = p, the diagonal term of 
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the corresponding matrix will vanish; i.e. the matrix coefficients of the matrix eigenvectors 
will vanish at p' = p. However there remain non-zero terms multiplying (f> nl {p ). These can 
be used as non-zero diagonal coefficient terms for the eigenvectors. The result is that one 
can obtain eigenvalues and eigenvectors directly from one’s matrix equation. As we shall 
see later this is no longer possible for the Linear potential. These techniques for the pure 
Coulomb case are also very well explained in references 13 and 14. 

However for the Linear potential this method does not work. (For the sake of simplicity 
let us discuss the / = 0 Linear potential equation only. The methods are identical for the 
higher l equation.) The reason that the above method does not work in Eq.(2.19) is 
because Q' 0 (y) has a double pole singularity and even after the subtraction, a principal 
value singularity is left. Thus the integral must be evaluated explicitly. However to do 
this we must know what the functions <j> n[) (p) are before we solve the problem ! The way 
around this dilemma is to expand <£„u (p) i n a suitable set of basis functions: 

M 

*„o (p) = £ CiJitp) (2.23) 


Inserting this expansion in Eq(2.19), multiplying by p 2 gj{p ) and integrating over p, we 
obtain: 


/ 1 *-gj{p)gi{p)dp+^r J J Q'M9j{p)(9i{p') - 9i{p)) d p' d p} 

* 0 0 0 

M 

= E n oY^Ci J p 2 9j{p)9i{p) d P ( 2 - 24 ) 


which is just the matrix equation: 


Y.AjiCi^E^Y.GjiCi (2-25) 

i i 

which is symmetric under interchange of t and j (equivalent to symmetry under interchange 
of p and p') thus ensuring that the eigenvalues are all real. 

The double integral still contains a principal value singularity. In order to treat this, 
the integral over p' is performed by integrating from 0 to 2 p and then 2p to oo with 
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the singularity at the midpoint of the first region, which is carried out using Gaussian 
quadrature with an even number of points. This type of integration yields the Cauchy 
principal value automatically 17 . When we solve (2.25) we get M eigenvalues E 10 to Emo 
and a corresponding set of M eigenvectors C x to Cm . Thus Eq.(2.25) is solved for the 
energies and the coefficients C;, which yield the wave function when substituted back 
into Eq.(2.23) Convergence is obtained by increasing the number of basis functions M and 
integration points. In order to obtain the wave function in coordinate space, one simply 
takes the Fourier transform < 7 j(r) of the basis functions <7,(p) and uses the same set of 
coefficients C,- but now multiplying </,-( r) to obtain the coordinate space wave function. 
(Thus it is very convenient to pick g,(p ) so that they have a simple Fourier transform.) 
For the masses and couplings considered in this paper a convenient set of functions p,(p) 
is 

< 7 ;(p ) = eip[— p 2 i 2 /M] (2.26) 

where M is the maximum number of functions used in the expansion Eq(2.23). Note 
however that for different masses and couplings 15 , a different set of basis functions is 
necessary to achieve rapid convergence. 

When solving the general Coulomb plus Linear problem one cannot take advantage 
of the simplicity of the Coulomb numerical procedure 13,14 by itself. One must employ the 
basis function expansion method described above. The basis functions appropriate to the 
Linear potential alone also turn out to be suitable for the general Linear plus Coulomb 
problem for the masses and couplings of this paper. 


2.4 NON-RELATIVISTIC RESULTS 


We have carried out many different tests of our methods. Firstly, for the pure Coulomb 
case we solved the problem with the method of references 13 and 14, which does not require 
any basis function expansion. We compared to the exact Coulomb energies and found that 
we could easily generate over 20 eigenvalues very accurately. Secondly, as an additional 
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check we also solved the pure Coulomb case using an appropriate set of basis expansion 
functions and were able to obtain about 10 eigenvalues quite accurately. Thirdly, the pure 
Linear problem was solved for l = 0 (see reference 15 for details) and compared to the exact 
results. (For the l = 0, pure Linear potential case, the exact eigenvalues can be obtained 
in terms of the roots of the Airy function). The calculated eigenfunctions also agreed with 
the exact results. Fourthly, the combined Coulomb plus Linear problem was solved with 
the expansion functions in Eq.(2.26) for / = 0, 1,2,3 and compared to a coordinate space 
calculation. (The coordinate space code integrates the Schrodinger equation out from the 
origin at r = 0 and in from large r, and matches the logarithmic derivatives at the classical 
turning point). Fifthly, the combined Coulomb plus Linear results were also compared to 
those listed in reference 18 and also with a coordinate space code. Excellent agreement 
was obtained. 

In summary, we have very thoroughly tested our methods for Coulomb plus Linear 
potentials for many partial waves against results from exact calculations, coordinate space 
codes and the results of other authors for both eigenvalues and wave functions. 


3. RELATIVISTIC TWO-BODY EQUATIONS AND INTERACTION KERNELS 


In traditional nuclear physics, the deuteron is the only two particle bound state system. 
It has been studied in both the nonrelativistic framework and also in numerous relativis- 
tic frameworks. Compared to the deuteron the qq system is a very rich system and its 
spectra provides an ideal testing ground in which a systematic study of the 3-dimensional 
relativistic equations can be made. 

The Bethe-Salpeter(BS) equation for the bound state problem in the center of mass 
frame is given by 


*(**) = (At / vW)C(^p 0 )*( P ,p 0 Ky (3.i) 

As mentioned above, there are infinitely many 3-dimensional reductions of the BS 
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equation. In this section we are going to work with six particular reductions which we 
believe to be a fair representative sample of the most commonly used 3-dimensional reduc- 
tions of the BS equation. In order to reduce Eq(3.1) in to a 3-dimensional equation, we 
replace the propagator G by a 3-dimensional propagator g which has the same elastic cut. 
A systematic study of these 3-dimensional relativistic equations for the problem of scatter- 
ing of scalar particles has been performed in reference 11. As stated in the introduction, 
in this paper we wiU make a similar study of the bound state of two particles interacting 
via a confining interaction. Some results have already been previously discussed 16 . The 
choice of the 3-dimensional propagator can be categorized into two types in general. One 
which renders the interaction to be instantaneous and one which does not. In this pa- 
per we study six 3-dimensional reductions, three of each type. Minimal Relativity (MR) 
equation 87 , Kadyshevsky (K) equation 9 , and Gross (G) equation 7 all of which retain the re- 
tardation in the interaction. The equations with instantaneous interaction(no retardation) 
are the Blankenbecler-Sugar (BBS)equation 8 , Kadyshevsky (KO) (without retardation) 9 
and Thompson (T) equation 10 . All six equations can be generically written as (compare 
to Eq.2.2) 


Di<f>(p) = - J dp'Vi{p,p')<j>{p') (3-2) 

where 0(p) is a Schrodinger like wave function. We will neglect the couplings to the 
negative energy channels since the subtraction method is the same for the coupled channel 
case. The D{ are given in table I and the index * can be MR, K, G, BBS, KO and T. Note 
that for MR, K and G equations the interaction V { has retardation and for the other three 
equations it does not. We will choose to use MR, BBS, K and KO equations to study the 
bound states of two scalar particles interacting via a confining interaction and G and T 

equations to study the bound state of spinor quarks . 

The confining interaction to be used in these relativistic equations is a straightforward 
generalization of the Linearly rising potential discussed in section 2. We simply replace 
the three vector q of Eq(2.2) by a four- vector q. Now q 2 is given by 

, ! =(p-p') 2 -(E p < 3 - 3 ) 
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In this generalization the form of the Coulomb type interaction and the confining interac- 
tion remain the same as in the nonrelativistic case but now q 2 is replaced by q 2 and the 
partial wave components of these interactions will be given by Eq(2.5) but for the equa- 
tions that include retardation (MR, K, G) the variable y is now replaced by y (instead of 
Eq.2.6 with rj = 0) where 


_ _ ( 3 . 4 ) 

2pp> 

Equations without retardation (BBS, K0, T) retain the original form of y in Eq.(2.6) with 
q = 0. Here p and p' are only the magnitude of the three vectors. Again we note that these 
relativistic interactions will introduce singularities as in the nonrelativistic case at q 0 
or at y = 1. Note also that although the variables are different the singularity structures 
are similar to the nonrelativistic case; i.e. the Coulomb interaction will have a logarithmic 
singularity and the confining interacting has higher order singularities. For the equations 
without retardation the interaction Vj'(Pi p ' ) is instantaneous and it is exactly the same as 
the nonrelativistic case. For the instantaneous interaction, relativistic effects come in to the 
equation only through the kinematics; i.e. only through the operator D{. The singularities 
in this interaction can be handled exactly the same way as in the nonrelativistic case. 

In the following subsections we will discuss how the singularities in the relativistic 
confining and Coulomb interactions can be treated properly. 


3.1 RELATIVISTIC COULOMB PROBLEM 


The relativistic generalization of the Coulomb interaction in the partial wave form is 
given by 


V, c (p',p) = — lim 

1 v 7 r 7—o pp' 

and by using the expression for Qi Eq(3.2) becomes 


(3.5) 
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oo oo 

D,Mp) + j-f ^ /<•!- MMpW = o (3.6) 

np o 0 

for the MR, K and G equations only. For the instantaneous equations BBS, KO, T instead 
of the above Eq.(3.6), we have simply the Schrodinger equation (2.14) but with the operator 
D{ replacing the Schrodinger propagator. Note that the only singularity in equation (3.6) 
arises from Q 0 (y). We want to handle this singularity in a similar fasion as in the nonrel- 
ativistic case; i.e. by adding and subtracting a term. But we must also be able to handle 
the added term analytically or numerically. Unfortunately because of the presense of re- 
tardation we cannot just subtract a 0 n j and use Eq.(2.12) as in the nonrelativistic case. In 
order to take advantage of Eq(2.12), we subtract a term propotional to the nonrelativistic 
interaction and obtain (compare to Eq. 2.15) 


Di<t>ni{p) + J [QM^nl{p')p' 

a 


Qo(y) 


pVnl(p) 

p'Pii y) 


Pi{y)dp 


+ — [p 2 <Mp)^t] - — [ w i-i{y) < t>ni{p)p'dp = o 

7T p l 7Cp J 

0 


(3.7) 


for the MR, K and G equations only. Again for the instantaneous equations BBS, KO, T 
instead of the above Eq.(3.7), we have simply the Schrodinger equation (2.15) but with 
the operator Z?,- replacing the Schrodinger propagator. Note that we again have at the 
singular point p' = p, y = 1 and Pj(y = 1) = 1 and by Taylor expanding Qu(y) around 
p' = p, one can show that the term in the square brackets vanishes at the singular point. 


3.2 RELATIVISTIC CONFINING PROBLEM 


In the case of the relativistic confining interaction, the functional structure of the in 
teraction is again the same as the nonrelativistic case but y replaced by y. We therefore use 
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the same type of subtraction used in the relativistic Coulomb case. That is, we subtract and 
add a term propotional to the nonrelativistic confining interaction. We obtain (compare 
to Eq. 2.22) 


OO 

d,Mp) + ^ J[Q'„iy)Mp') - 


P':(y)dp' 

0 

o° 

- ^ = 0 (3 ' 8) 

o 

for the MR, K and G equations only. Once more for the instantaneous equations BBS, 
KO, T instead of the above Eq.(3.8), we have the Schrodinger equation (2.22) but with the 
operator D{ replacing the Schrodinger propagator. The factor E 2 /m 2 in the subtracted 
term of the first integral in Eq(3.8) is necessary in order to cancel the singularity arising 
from Q',j(y) exactly at the singular point. This can easily be seen by Taylor expanding 
Q'u(y) at p' = p. Equation (3.8) is now ready to be solved for various choices of D{ when 
there is retardation in the interaction. For cases without retardation Vi is identical to 
the nonrelativistic problem and the subtraction technique developed in the nonrelativistic 
section can be used. 


3.3 RELATIVISTIC RESULTS AND CONCLUSIONS 


The main purpose of the present paper is to present the theoretical subtraction tech- 
niques necessary to solve two-body relativistic bound state equations in momentum space. 
Therefore equations (3.7) and (3.8) are our major results. 

Nevertheless for the sake of illustration we shall present some numerical solutions for 
the pure confining problem with equation (3.8) written in terms of a single channel. Such 
results will at least allow us to see whether our theoretical methods give reasonable results. 
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The usefulness of these relativistic equations depends on the extent to which they reproduce 
global properties of the spectrum characterized by the dependence of the energy E n i on 
the principal quantum number n. This dependence is most easily revealed by studying the 
ratio E n i/Eu . E n i is related to the total energy W n i through E n \ = W n i — 2m. Tables 
2, 3 and 4 contain the results for the ratio E n i/En for the equations listed above for a 
reasonable choice of mass and coupling parameters. I values range from 0 to 2. 

There are three observations to make from these tables. First all of the energy ratios 
are reasonably close to the non-relativistic results for heavy quark masses. Second the 
difference between the relativistic results and the non-relativistic results gets bigger for 
smaller quark mass. Third, the higher radial excitations show more pronounced relativistic 
corrections, which is consistent with the virial theorem 3 for a positive power law potential 
which requires larger kinetic energies for orbits with greater average radii. These results 
lead us to conclude that our theoretical methods are valid and give us confidence that 
the methods developed herein will be suitable when a full coupled channels calculation is 
performed and compared to experimental data. 

In conclusion we have presented the theoretical subtraction techniques necessary to 
solve two-body relativistic bound state equations in momentum space with Coulomb plus 
Confining interactions. Future work will be devoted to including spinors and coupling to 
the negative energy channels in all six equations so that detailed comparisons to experiment 
can be carried out. 
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Table 1 

Di operators for relativistic equations 
G and T equations are describing pseudoscalar mesons with spinor quarks 
The other four relativistic equations are for scalar quarks 


i 

Name 

Di 

Retardation 

MR 

Minimal 

Relativity 

4 E k {E k 2 - W 2 / 4) 

Yes 

BBS 

Blanckenbecler 

Sugar 

same as MR 

No 

K 

Kadyshevsky 

2 E k 2 {E k - W/2) 

Yes 

K0 

Kadyshevsky 

same as K 

No 

G 

Gross 

2 E k - W 

Yes 

T 

Thompson 

same as G 

No 
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Table 2 


Energy ratios for pure Confining interaction with l = 0. 

G and T equations are for spinor quarks with k = 0.2GeV 2 . The other four relativistic 
equations are for scalar quarks with k = O^GeV^. The nonrelativistic(NR) equation is 
with k = 0.2GeF 2 . All masses are in units of GeV. 


n 

MR 

BBS 

K 

K0 

G 

T 

NR 

Mass 

1 

1.73 

1.71 

1.74 

1.72 

1.79 

1.72 

1.75 

1.5 

2 

2.31 

2.27 

2.34 

2.30 

2.47 

2.30 

2.36 

1.5 

3 

2.81 

2.75 

2.87 

2.80 

3.09 

2.80 

2.90 

1.5 

1 

1.58 

1.50 

1.68 

1.54 

1.90 

1.67 

1.75 

0.5 

2 

2.00 

1.82 

2.21 

1.89 

2.73 

2.18 

2.36 

0.5 

3 

2.35 

2.08 

2.65 

2.16 

3.52 

2.62 

2.90 

0.5 

1 

1.51 

1.41 

1.66 

1.44 

1.98 

1.63 

1.75 

0.3 

2 

1.87 

1.65 

2.13 

1.69 

2.92 

2.11 

2.36 

0.3 

3 

2.17 

1.84 

2.52 

1.89 

3.84 

2.51 

2.90 

0.3 
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Table 3 


ri 

Energy ratios — for pure Confining interaction with 1 = 1 . 
Notation and units are the same as Table 2. 


n 

MR 

BBS 

1 

1.44 

1.43 

2 

1.82 

1.79 

3 

2.16 

2.11 

1 

1.38 

1.31 

2 

1.67 

1.54 

3 

1.92 

1.72 

1 

1.37 

1.27 

2 

1.64 

1.45 

3 

1.87 

1.60 


K 

K0 

G 

1.45 

1.44 

1.49 

1.85 

1.81 

1.92 

2.20 

2.15 

2.33 

1.47 

1.35 

1.56 

1.84 

1.59 

2.09 

2.17 

1.79 

2.60 

1.52 

1.30 

1.61 

1.90 

1.50 

2.21 

2.22 

1.65 

2.80 


T 

NR 

Mass 

1.43 

1.45 

1.5 

1.80 

1.85 

1.5 

2.13 

2.20 

1.5 

1.39 

1.45 

0.5 

1.71 

1.85 

0.5 

2.00 

2.20 

0.5 

1.36 

1.45 

0.3 

1.67 

1.85 

0.3 

1.93 

2.20 

0.3 
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Table 4 


Energy ratios for pure Confining interaction with 1 = 2 . 

Notation and units are the same as Table 2. 


n 

MR 

BBS 

K 

K0 

G 

T 

NR 

Mass 

1 

1.32 

1.31 

1.33 

1.31 

1.35 

1.31 

1.33 

1.5 

2 

1.59 

1.57 

1.62 

1.59 

1.68 

1.58 

1.62 

1.5 

3 

1.84 

1.81 

1.88 

1.84 

1.99 

1.82 

1.89 

1.5 

1 

1.29 

1.23 

1.37 

1.26 

1.41 

1.27 

1.33 

0.5 

2 

1.52 

1.40 

1.67 

1.45 

1.80 

1.51 

1.62 

0.5 

3 

1.71 

1.55 

1.94 

1.61 

2.19 

1.72 

1.89 

0.5 

1 

1.30 

1.20 

1.44 

1.23 

1.45 

1.25 

1.33 

0.3 

2 

1.52 

1.35 

1.78 

1.39 

1.89 

1.47 

1.62 

0.3 

3 

1.72 

1.47 

2.06 

1.52 

2.33 

1.67 

1.89 

0.3 
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